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Abstract: This note documents the use of the POWHEG-hvq package, a generator for 
heavy flavour hadroproduction at next-to-leading order in QCD, that can be easily inter- 
faced to shower Monte Carlo programs, in such a way that NLO and shower accuracy are 
both maintained. 
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1. Introduction 

The POWHEG-hvq program is a liard event generator for heavy quark production in 
hadronic collisions. It is accurate at the next-to-leading order in QCD (NLO from now 
on), and it can be interfaced to Shower Monte Carlo programs (SMC from now on) like 
HERWIG and PYTHIA, in such a way that both the leading logarithmic accuracy of the 
shower and the NLO accuracy are maintained in the output. It is thus an alternative to 
the MC@NLO heavy flavour production program of ref. [1]. The code can be found in 

http : //moby . mib . inf n . it/~nasoii/POWHEG. 

The program is an implementation of the heavy flavour NLO cross sections of refs. [2, 3, 4], 
according to the formalism of refs. [5] and [6]. A detailed description of the implementation 
is given in ref. [7]. In the case of tt production, spin correlations in top decays are included 
with a method similar to the one discussed in ref. [8] . The relevant matrix elements for the 
spin correlations are the same used in the MC@NLO package [9], and were obtained using 
MadEvent [10]. Our decay package takes into account the finite width of the top quarks 
and of the W. 

In this note, we give all the necessary information to run the program. 
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2. Installation 

The program comes as a tarred-gzipped file POWHEG-hvq.tar .gz. It can be installed with 
the following commands 

$ tar -zxf POWHEG-hvq.tar.gz 
$ cd POWHEG-hvq 
$ make <target> 

where the choice of the target depends upon the way one wants to interface the pro- 
gram with a Shower Monte Carlo implementation. The Makefile is set up to use the 
compiler g77 on Linux platforms. If one wishes to use gf ortran, one should uncomment 
the appropriate lines in the Makefile. Notice that some versions of gf ortran may not 
support the idate and time intrinsics. These are used in the mbook.f file, that in turn 
is used in the examples. Thus, one may also need to comment out the calls to idate and 
time in order to run the examples with gf ortran. 

3. Modes of operation 

The program POWHEG-hvq generates hard events. The hard events can then be fed 
into an SMC program for subsequent showering. POWHEG-hvq saves the hard event 
information according to the Les Houchcs Interface for User Processes (LHIUP from now 
on) conventions [11]. The SMC should also comply with these conventions (as is the case 
for PYTHIA and HERWIG) in order to be used in conjunction with POWHEG-hvq. 
The program can be run in three ways. 

• POWHEG-hvq generates hard events, and store them in a file. An SMC reads the 
file and showers it. 

• POWHEG-hvq is linked directly to the SMC. In this case the events are generated 
and immediately showered, without intermediate storage. 

• POWHEG-hvq is run as a standalone program, and the produced hard events are 
analyzed without showering. The output yields in this case NLO distribution with 
LL resummation of soft gluon effects. 

3.1 Storing the user events 

The easiest way to interface POWHEG-hvq to an SMC is to simply store the hard events 
in a file (which we call the event file) , and in a subsequent run read the events and pro- 
cess them with the SMC. The format of the event file supported by POWHEG-hvq is the 
"Standard format for Les Houches event files", documented in ref. [12]. The program for 
the generation of the Les Houches Event Files (LHEF from now on) can be built with the 
command 
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$ make main-lhef 

The event file is named pwglhef .events (the user is given the possibihty to change the 
file name, as documented in the next section). An example program that reads the event 
file, showers it with HERWIG and analyzes it the can be built as follows 

$ make main-HERWIG-lhef . 

A similar program, named main-PYTHIA-lhef is provided for PYTHIA. The HERWIG or 
PYTHIA program should be installed in the POWHEG-hvq directory. In the case of HER- 
WIG, the appropriate include files should also be present. As can be evinced from the Make- 
file, the fortran files relevant for these examples are main-HERWIG-lhef . f, herwig6510.f 
(main-PYTHIA-lhef .f, pythia6326.f for PYTHIA), mbook.f and analize-hvq. f . The 
file analize-hvq. f is a minimal analysis program, provided as a starting example for more 
complex analysis. It is adequate for both HERWIG and PYTHIA (since it uses the stan- 
dard common blocks of ref. [13]). It uses the histogramming package of M.L. Mangano, 
mbook.f, and it produces topdrawcr output in the file pwgoutput .top. It can be used for 
charmed and bottomed hadron production, and for tt production. For charm and bottom, 
it analyses a few kinematic obscrvablcs, looking only at D^, D^, D for charm, and B^,B^, 
B for bottom. For top pair production, the program looks for lepton pairs with transverse 
momenta above a given cut. 

POWHEG-hvq also provides an output that is as close as possible to the MC@NLO 
user file format, so that users that have already implemented MC@NLO in their analysis 
frameworks should be able to run POWHEG-hvq interfaced to HERWIG with no extra 
work. The corresponding Makefile target is main-mcatnlof 1. The event file name is 
pwgmcatnlof mt . events. 

3.2 Interfacing POWHEG-hvq with a shower Monte Carlo program 

One should create a main program that initializes the SMC to make it ready to accept a 
user process, and provide the following routines 

subroutine UPINIT 
call pwhginit 
end 

subroutine UPEVNT 
call pwhgevnt 
end 

that are the only link to the POWHEG-hvq program. The main program should call 
the appropriate subroutines to run the SMC. If the SMC is compliant with the LHIUP, it 
will call the routines UPINIT and UPEVNT in order to initialize, and to generate the hard 
events. The routine pwhinit performs the initialization of POWHEG-hvq, setting up all 
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the grids that are necessary for the efficient generation of the events, and it also initiahzes 
the process common block of the LHIUP. Each call to pwhgevnt results in the generation 
of one event, and its storage in the LHIUP event common block. 

When using HERWIG, one must remove the dummy subroutines UPINIT and UPEVNT 
that are present in the HERWIG source file. 

We provide two examples, that can be built with the commands 

$ make main-HERWIG 
and 

$ make main-PYTHIA 

The analysis program is the same one as that described in section 3.1. 

3.3 POWHEG-hvq as a standalone program 

In this case, the main program should have the structure 

program MAIN 
call pwhginit 
do j=l,NEVENTS 

call pwhgevnt 
c call some analysis routines here 

enddo 

c Print out results 
end 

No examples are provided. The analysis routines should make use of the information 
stored in the LHIUP, as documented in ref. [11]. 

4. Input parameters 

POWHEG-hvq provides an independent facility to set the input parameters for the run. 
All parameters are stored in a file, named powheg. input. The format of the file is as 
follows 

1. Lines are no more than 100 characters long. 

2. Empty (blank) lines are ignored 

3. If a # or a ! appears at any point in a line, the part of the line starting from the # or ! 
symbol up to its end is blanked. 
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4. An entry has the format: 
name value 

usually followed by a ! and a comment to clarify the meaning of the variable. The 
name keyword has no more than 20 characters, and value is an integer or floating point 
number. 

5. A maximum of 100 keywords are allowed. 

If the file powheg . input is not present, the program asks the user to enter a prefix, and then 
looks for the file <prefix>-powheg. input. In this case, all the files created by POWHEG- 
hvq in the current run will carry the prefix <pref ix>- instead of pwg. 

The input parameters are read by the (real * 8) function powheginput (string) , in 
file powheginput .f . The statement 

rvalue=powheginput ( 'myparm' ) 

returns the value of token myparm stored in powheg. input. If the token is not found in the 

input file, a message is printed, and the program is stopped. The file is read only once, on 
the first invocation of the function powheginput, and token- value pairs are stored in inter- 
nal arrays, so that subsequent calls to powheginput are relatively fast. With the statement 

rvalue=powheginput ( ' #myparm ' ) 

in case the token myparm is not present, the program does not stop, and returns the 
value —10^. The file powheginput . f is a standalone code, and can be linked to any 
program. In this way, an SMC that is reading an event file may get parameters of the 
POWHEG-hvq run, if it needs too. 

We document here a typical input file powheginput.dat: 

! Heavy flavour production parameters 

maxev 500000 ! number of events to be generated 
ihl 1 ! hadron 1 type 
ih2 -1 ! hadron 2 type 
ndnsl 191 ! pdf for hadron 1 
ndns2 191 ! pdf for hadron 2 

The integer ihl,ih2 and ndnsl, ndns2 characterize the hadron type and PDF used in 
POWHEG-hvq. The numbering scheme is documented in the file hvqpdfpho.f. In that 
file, in the routine PRNTSF, all PDF sets available and the corresponding set number are 
listed. In particular, 191 corresponds to MRST2002. The hadron type in ihl,ih2 is for a 
nucleon (i.e. the average of a proton and a neutron), 1 (-1) for a proton (antiproton), 2 (-2) 
for a neutron (antincutron), and 3 (-3) for a tt"*" (vr~). Thus ihl=l,ih2=-l corresponds to 
proton-antiproton collisions. 
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ebeaml 980 ! energy of beam 1 
ebeaiii2 980 ! energy of beam 2 



We assume that beam 1 and 2 move along the third axis in the positive and negative 
direction respectively. 

qmass 175 ! mass of heavy quark in GeV 

facscfact 1 ! factorization scale factor: muf act=muref *f acscf act 
renscfact 1 ! renormalization scale factor: muren=muref *renscf act 
bbscalevar 1 ! if use muref=qmass in Bbar calculation 

Factorization and renormalization scale factors appearing here have to do with the com- 
putation of the inclusive cross section (i.e. the B function [5, 6, 7]), and can be varied by 
a factor of order 1 in both directions to study scale dependence. Normally the reference 
scale is set equal to the transverse mass of the heavy quark in the rest frame of the qq 
system. If bbscalevar is set to 0, the reference scale is chosen equal to the heavy quark 
mass. Other choices require a modification of the the subroutine setscalesbb in the file 
phy spar-hvq . f . 

The following is only needed if the quark is a top, and we want it to decay with the inclu- 
sion of spin correlations 



topdecaymode 20000 ! an integer of 5 digits representing the decay mode. 

The value of the token is formed by five digits, each representing the maximum num- 
ber of the following particles in the (parton level) decay of the tt pair: e^, fi^, r^, u^, c^. 
Thus, for example, 20000 means the t e'^Ugb, t e~vj), 22222 means all decays, 10011 
means one of the top quarks goes into electron or antielectron, and the other goes into 
any hadron, 00022 means fully hadronic, 00011 means fully hadronic with a single charm, 
00012 fully hadronic with at least one charm. If all digits are 0, neither the t nor the i 
are decayed. Values that imply only one t decay (for example 10000) are not implemented 
consistently. 

In case topdecaymode is different from 0, more parameters are needed for the decay kine- 
matics, and are used exclusively for decays 



tdec/wmass 80.4 
tdec/wwidth 2.141 
tdec/bmass 5 
tdec/twidth 1.31 



W mass for top decay 
W width 

b quark mass in t decay 
top width 



tdec/elbranching 0.108 ! W electronic branching fraction 
tdec/emass 0.00051 ! electron mass 
tdec/mumass 0.1057 ! mu mass 
tdec/taumass 1.777 ! tau mass 



-6- 



tdec/dmass 0.100 
tdec/umass 0.100 
tdec/smass 0.200 
tdec/cmass 1.5 



u mass 



d mass 



s mass 



charm mass 



tdec/sin2cabibbo 0.051 ! sine of Cabibbo angle squared 

Spin correlations in the decay are implemented, and effects due to the finite width of 
the top and of the W are also accounted for. 

The following parameters control the operation of the POWHEG-hvq program: 

! Parameters to allow-disallow use of stored data 
use-old-grid 1 
use-old-ubound 1 

The meaning of these tokens requires a little knowledge of the operation of POWHEG- 
hvq. Before the program starts generating events, the integral of the inclusive cross section 
is computed, and a grid is set up for the generation of Born-like configurations. Similarly, 
in the generation of hard radiation a grid is computed to get an upper bounding function 
to the radiation probability. The generation of the grids is time consuming, but the time 
spent in this way is negligible in a normal run with hundreds of thousands of events be- 
ing generated. On the other hand, it is useful (for example, when debugging an analysis 
program) to skip the generation stage. For this purpose, the grid for the generation of 
Born-like kinematics is stored in the file pwggrid.dat. 

If use-old-grid is set equal to 0, and pwggrid.dat exists and is consistent, it is loaded, 
and the old grid and old value of the cross section are used. Otherwise, a new grid is 
generated. Observe that the program does check the file for consistency with the current 
run, but the check is not exhaustive. The user should make sure that a consistent grid is 
used. 

The token use-old-ubound has the same role as use-old-grid, but it applies to the upper 

bounding array that is used in the generation of radiation. 

The following parameters are used to control the grids generation 

! parameters that control the grid for Born variables generation 

ncalll 10000 ! number of calls for initializing the integration grid 

itmxl 5 ! number of iterations for initializing the integration grid 

ncall2 100000 ! number of calls for computing integral and upper bound grids 

itmx2 5 ! number of iterations for computing integral and upper bound grids 

foldx 1 ! number of folds on x integration 

foldy 1 ! number of folds on y integration 

foldphi 1 ! number of folds on phi integration 

! Parameters that controll the generation of radiation 

nubound 100000 ! number of bbarra calls to setup upper bounds for radiation 
iymax 1 !<=10, number of intervals in y grid to compute upper bounds 
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ixmax 1 !<=10, number of intervals in x grid 

xupbound 2 ! increase upper bound for radiation generation by given factor 
The values of some of the tokens may be changed in the following cases: 

• If the integration results have large errors, one may try to increase ncalll, itmxl, 
ncall2, itiiix2. 

• If the fraction of negative weights is large, one may increase f oldx, f oldy, f oldphi. 
Allowed values are 1, 2, 5, 10, 25, 50. The speed of the program is inversely propor- 
tional to the product of these numbers, so that a reasonable compromise should be 
found. 

• If there are too many upper bound violation in the generation of radiation (see Section 
6), one may increase nubound, and/or xupbound. 

• If the efficiency in the generation of radiation is too small, one may try to increase 
iymax, ixmax. 

In oder to check whether any of these conditions occurs, the user should inspect the file 
pwgstat . dat at the end of the run, as illustrated in sec. 6. 

5. Examples 

Examples of powheg. input files are given in the directories c-tev, b-tev, t-tev, tdec-tev 
and c-lhc, b-lhc, t-lhc, tdec-lhc. In the examples in the tdec-tev, tdec-lhc di- 
rectories the ti pair is decayed semileptonically into an electron and an antielectron by 
POWHEG-hvq. In the t-tev, t-lhc directories the top decay is handled by the SMC, 
according to its own parameters, but without including spin correlations. In all examples, 
the choice of the parameters that control the grid generation is such that a reasonably small 
fraction of negative weights is generated. Even in the extreme case of charm production at 
the LHC, where a fraction of negative weight less than 3% is achieved. 

6. Counters and statistics 

Several results relevant to the interpretation of the output of the run are written to the 
file pwgstat.dat. The fraction of negative weights, the total cross section, the number of 
upper bound failures in the generation of the inclusive cross section, and the generation 
efficiency, together with failures and efficiency in the generation of hard radiation, are 
printed there. These numbers are sufficient to take action in case of problems. 

A call to the subroutine printstat causes a printout of all POWHEG-hvq counters in 
pwgstat . dat file. 
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7. Using the PDF sets 



POWHEG-hvq uses the PDF implementation of hvqpdfpho.f . In that file, in the routine 
PRNTSF all pdf sets available, and the corresponding set numbers are listed. In order to 
use a set, the corresponding data file must be present in the directory where the run is 
performed. If the file is not found a message is printed with the name of the missing file, 
and the program stops. In this case one should copy the missing file from the directory 
pdf data to the current directory, or set up a symbolic link to it. 

8. Random number generator 

POWHEG-hvq uses the RM48 random number generator, documented in the CERNLIB 
writeups. This generator has default initialization. If a user wishes to start the program 
with different seeds, he should call rm48in(iseed,nl ,n2) (before calling the pwhginit rou- 
tine) in order to seed the generator with the integer iseed, and skip the first nl+n2*10**8 
numbers, as documented in the CERNLIB manual. This affects the POWHEG-hvq ran- 
dom number sequence. If the program is interfaced to an SMC, the user should also take 
care to initialize the seeds of the latter. 
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